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摘 要 :冰川 径流 是 西北 干旱 区 径流 的 主要 组 成 部 分 ,研究 未 来 气候 变化 对 冰川 径流 的 影响 对 西北 
干旱 区 径流 至 关 重 要 。 以 博 尔 塔 拉 河 上 游 源 流 区 为 研究 区 ,构建 铭 入 冰川 模块 的 SWAT 模 型 ,模拟 
温泉 水 文 站 1972 一 2018 年 月 径流 过 程 ,并 在 此 基础 上 研究 了 气候 变化 情景 下 (RCP4.5 和 RCP8.5) 
未 来 (2020 一 2050 年 ) 气 候 变 化 对 冰川 径流 的 影响 。 结 果 表 明 :SWAT 模 型 能 够 很 好 地 模拟 源流 区 
径流 变化 过 程 ,在 整个 模拟 期 间 , 径 流 数据 的 纳什 系数 (NSE) 为 0.82, 偏 差 百分比 (PBIAS ) 为 -=3.22% , 
均 方 根 误差 与 实测 值 标准 差 的 比值 (RSR ) 为 0.42 ,决定 性 系数 ( 尼 ) 为 0.84, 模 型 性 能 评定 为 优 。 根 
42 CMIPS 气候 模式 2 种 情景 的 模拟 结果 ,2 种 情景 模拟 未 来 总 径流 都 呈现 出 增加 趋势 ,分 别 将 以 
0.31x10 m°: (10a) 和 0.40x10m (10a) 的 速度 继续 增加 ,冰川 径流 占 比 较 历 史 时 期 的 27.619% 分 
别提 升 了 4.84% 和 9.38%。 冰 川 径流 增加 是 径流 量 增加 的 主要 原因 。 通 过 相关 性 分 析 发 现 , 随 着 
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气温 的 升 高 ,冰川 消融 时 间 提 前 ,冰川 消融 加 速 ,冰川 积累 时 间 减 少 ,导致 冰川 面积 进一步 的 缩 
减 。 研 究 结果 可 为 博 河 地 区 水 文 资料 历史 变化 、 未 来 注 变 趋势 和 预防 气候 变化 带 来 的 潜在 风险 提 


供 依 据 。 


X 键 i: SWAT 模 型 ;气候 变化 ; CMIP5 气候 情景 模式 ; 冰川 径流 ; 博 尔 塔 拉 河 


文章 编号 : 


艾 比 湖 流域 (79°53'~83°53'E,44°02'~45°23'N) 
位 于 新 疆 西 北部 ,是 准 踢 尔 盆 地 西部 最 低 的 汇 水 中 
心 。 自 20 世 纪 50 年 代 以 来 ,由 于 大 规模 的 水 土 开 
发 ,灌区 人 口 灌溉 面积 和 引水 量 大 幅度 增加 ,原来 
有 水 注入 艾 比 湖 的 大 部 分 河流 先后 断 流 ,导致 人 湖 


1000 - 6060(2023)09 - 1432 - 11(1432 ~ 1442) 


编目 条 数 已 经 退 减 到 243 条 ,面积 消融 至 113.75 kim’. 
两 次 冰川 编目 数据 相 比 , 博 河 冰川 消失 12 条 ,面积 
缩减 103.15 km 。 随 着 全 球 进一步 增 暖 ,气候 模式 
预 估 全 球 水 循环 进一步 增强 ,水 循环 季节 差异 增 大 … 。 
博 河上 游 地 区 位 于 高 山区 ,在 气候 变化 背景 下 , 山 


水 量 急 剧 减 少 ,流域 径流 不 仅 要 承担 人 们 生产 生活 
供水 任务 ,还 要 为 维持 当地 良好 生境 提供 生态 用 
水 ,2019 年 生态 环境 用 水 量 仅 占 总 用 水 量 的 3.1%， 
当地 社会 经 济 - 生 态 环境 用 水 矛盾 日 益 增加 。 博 尔 
塔 拉 河 (简称 博 河 ) 是 艾 比 湖 的 主要 入 湖水 源 , 近 40 a 
补给 水 量 达 入 湖 总 量 的 70%, 博 河水 量变 化 直接 关 
平 艾 比 湖面 积 增 减 及 湖 周 生 态 状况 ,对 艾 比 湖 生态 
修复 保护 发 挥 基 础 性 的 支撑 作用 。 博 河 第 一 次 冰 
川 编目 条 数 为 255 条 、 面 积 216.91 km ;第 二 次 冰川 


收 稿 日 期 : 2022-10-24; ft] ABA: 2022-12-14 


区 冰川 消融 加 剧 ,导致 极端 水 文 事件 的 发 生 概率 增 
多 ,河流 水 文 过 程 将 会 变 得 更 为 复杂 ”“。 因 此 ,开展 
气候 变化 背景 下 博 河 径流 变化 预测 ,对 预防 气候 变 
化 带 来 的 潜在 水 资源 风险 ,保障 当地 社会 经 济 发 展 
与 生态 文明 建设 显得 十 分 迫切 。 

目前 ,针对 径流 的 研究 主要 包括 寻找 数据 间 的 
最 优 数 学 关系 为 目标 的 数据 驱动 模型 和 以 水 文学 
概念 为 基础 考虑 物理 过 程 的 水 文 模型 "。 与 数据 驱 
动 模型 相 比 ,水 文 模型 能 真实 有效 地 模拟 流域 径 
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流 形成 过 程 ,客观 地 反映 气候 或 下 垫 面 变化 对 流域 
水 循环 的 影响 程度 ,是 研究 水 文 循 环 机 理 和 过 程 、 
解决 复杂 水 文 问题 的 有 效 基础 工具 “。 在 众多 水 文 
模型 中 ,SWAT 模 型 具有 很 高 的 普 适 性 ,在 各 流域 的 
模拟 应 用 中 ”7 , 均 取 得 了 较 好 的 效果 ;在 对 气候 变 
化 的 研究 过 程 中 ,SWAT 模 型 成 为 国内 外 研究 人 员 
开展 工作 最 常用 的 模型 之 一 ,如 Shukla 等 2 Cao | 
Zhao 等 中 分 别 运 用 SWAT 模 型 在 Satluj 流域 . 呼 伦 
湖 流域 .白山 使 地 等 开展 了 气候 变化 情景 对 径流 的 
影响 研究 。 同 时 ,SWAT 模 型 和 其 他 模块 具有 耦合 
性 , 杨 明智 等 5 基于 SWAT 模 型 柑 入 水 资源 配置 模 
块 ,实现 了 “自然 -社会 ”二 元 水 循环 与 水 资源 实时 
配置 一 体 化 模拟 功能 ; Yin 等 将 SWAT 模 型 嵌入 冰 
川 模 块 ,解决 了 SWAT 模 型 在 高 寒山 区 模拟 精度 不 
高 的 问题 ,并 定量 估算 了 叶 尔 羌 河 流域 冰川 对 径流 
的 贡献 。 可 以 看 出 ,SWAT 模 型 不 仅 在 广泛 应 用 于 
各 流域 和 气候 变化 的 研究 中 ,还 能 针对 自身 薄弱 
点 ,因地制宜 地 耦合 相关 模块 ,提高 模型 模拟 精度 。 

在 博 河 以 往 的 研究 中 ,多 侧重 从 融雪 描述 径流 
机 理 ” AK a) aK EO) KEk. .地 下 地 表 水 
同位 素 变化 特征 等 方面 的 研究 ,很 少 关注 气候 变 
化 对 博 河 流域 水 文 过 程 的 影响 。 本 研究 以 博 河 上 
游 为 研究 区 ,将 冰川 模块 朋 人 SWAT 模型 中 ,构建 能 
够 刻画 冰川 径流 的 冰川 增强 型 SWAT 模 型 ,并 耦合 
CMIP5 气 候 变化 模式 中 的 不 同情 景 ,预测 未 来 径流 
增 减 情况 ,揭示 引起 径流 变化 的 主要 驱动 因子 ,以 
期 为 变化 环境 下 的 流域 水 资源 管理 措施 提供 辅助 
支撑 。 
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1 研究 区 概况 


博 河 流域 (79°53'~83°53 禄 ,44°02'~45°23'N) 地 
处 西天 山北 坡 的 阿拉 套 山 ,全 长 252 km, 近 60 a F 
均 径 流量 为 $.1x108m ,经 温泉 县 和 博 乐 市 后 注入 艾 
比 湖 ,并 以 温 录 县 和 博 乐 市 为 界 ,将 博 河 分 为 上 、 
中 下 游 , 其 中 冰川 面积 为 110.3 km ,冰川 补给 量 约 
为 1.05x10m (图 1)。 博 河 温泉 站 以 上 集 水 面积 ， 
2206 km ,全 长 97 km。 年 均 气 温 为 6.6~7.8 % ,其 中 
5 一 9 月 的 气温 达到 10 % 以 上 ,最 高 是 7 月 ;多 年 平 
均 降 水 量 为 288.63 mm, 主要 集中 在 4 一 9 月 ,6 个 月 
内 降水 量 约 占 全 年 的 80.7% ,其 中 降水 量 最 大 的 月 
份 是 7 月 。 可 见 ,夏季 是 气温 增高 降水 增多 的 集中 
时 段 , 最 高 月 份 均 出 现在 7 月 , 雨 热 同 期 。 全 年 径流 
呈现 中 间 高 两 头 低 的 单 峰 分 布 ,年 内 径流 主要 集中 
发 生 在 6 一 8 月 ,3 个 月 内 径流 量 约 占 全 年 的 40.62% , 
其 他 9 个 月 份 径流 量 相 对 平均 , 均 为 0.2x10: m? 
左右 。 


2 数据 与 方法 


2.1 研究 方法 

2.1.1 SWAT 模 型 SWAT 模 型 是 在 美国 农业 部 
(USDA) 的 农业 研究 中 心 Jeff Arnold 博士 的 主持 下 ， 
于 1994 年 将 SWRRB 模型 与 ROTO 模型 整合 到 一 
起 ,开发 出 了 SWAT 模 型 。 该 模型 经 过 几 十 年 的 发 
展 ,已 经 被 广泛 应 用 于 径流 模拟 VED ,点 源 污染 等 
方面 下 。 在 SWAT 模 拟 各 计算 单元 的 水 循环 过 程 


博 尔 塔 拉 河 
人 | 湖泊 


图 1 艾 比 湖水 系 分 布 


Fig. 1 Distribution of Ebinur Lake river System 
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时 ,水 量 平衡 法 则 是 根本 驱动 力 。 模 型 首先 根据 数 
字 高 程 模型 (DEM ) 数 据 、 土 壤 类 型 .土地 利用 类 型 
及 坡度 , 细 分 流域 反映 不 同 作物 和 土壤 类 型 的 燕 散 
发 差异 ,然后 模拟 各 水 文 响应 单元 的 径流 量 , 最 后 
汇总 得 到 流域 的 总 径流 量 ,物理 意义 明确 。 


SW, =SW, + $ (Rag -Qua -E -Wup Oye) (1) 
isi 


式 中 : SW, 为 土壤 最 终 含水 量 (mm); SW, 为 第 i 天 
的 土壤 初始 含水 量 (mm);/ 为 时 间 (d); Ry, 为 第 ;天 
的 降水 量 (mm); 0, 为 第 i 天 的 地 表 径流 量 (mm); 
E, 为 第 ;天 的 蒸 散发 量 mm); 到 .为 第 ;天 从 土壤 剖面 
进入 包 气 带 的 水 量 Amm); 0, 为 第 i 天 回归 流 的 水 量 
(mm). 
2.1.2 冰川 模块 算法 SWAT 模 型 具有 较 强 的 径流 
演化 与 溶质 迁移 模拟 功能 ,但 对 干旱 区 内 陆 河 流域 
径流 的 重要 组 分 一 一 冰川 融 水 模拟 考虑 不 足 。 
Yin 等 “根据 冰川 融化 因子 和 气温 的 关系 ,考虑 太 
阳 辐 射 和 地 形 等 因素 的 影响 ,修正 了 温度 指数 模 
型 ,对 冰川 水 文 过 程 的 模拟 更 具 物 理 意 义 。 
M- | ut Riel ia 人 三 了 ae 了 > T iiie 
0 jp LST rie 
式 中 :为 冰川 日 融化 量 (mm); 为 冰 温 融化 因 
F; Re 为 冰 辐 射 融 化 因子 ; oa 为 潜在 大 阳 辐 射 ;7 
为 日 平均 气温 (% ); Ti 为 冰川 融化 温度 闵 值 
(C)。 有 关 应 用 表明 ”” ,这 一 改进 对 提高 干旱 内 
陆 河 径流 模拟 具有 较 好 的 完善 作用 。 
2.1.3 模型 评价 本 文选 定 4 个 评价 指标 作为 评判 
模拟 结果 好 坏 的 标准 ,分 别 是 纳什 系数 (NSE) 、 均 方 
根 误 差 与 实测 值 标 准 差 的 比值 (RSR) 偏差 百分比 
(PBIAS) 和 决定 性 系数 (RR) ,其 中 尼 表 示 模 拟 值 与 
实测 值 的 相关 程度 , 取 值 范 围 为 0~1, 越 接近 1 表明 
相关 性 越 好 ,具体 评价 标准 ” 见 表 1。 各 指标 计算 
公式 如 下 : 


(2) 


(3) 


Fel 模型 评价 标准 


Tab.1 Model evaluation criteria 


模拟 结果 NSE RSR PBIAS/% 
优 0.75<NSE<1.00 0.00<RSR<0.50 IPBIASI<10 

良 0.65<NSE<0.75 0.50<RSR<0.60 10<IPBIASI<15 
合格 0.50<NSE<0.65 0.60<RSR<0.70 15<IPBIASI<25 
不 合格 NSE<0.50 RSR>0.70 IPBIASI=25 


HE: NSE 为 纳什 系数 ;RSR 为 均 方 根 误差 与 实测 值 标准 差 的 比值 ; 
PBIAS 为 偏差 百分比 。 下 同 。 


[Ser =o) 0% cg 


Rs (6) 


n 


2o" = om) Do” - in) 

式 中 : Qe 为 时 间 步 长 ;内 模型 模拟 的 河道 径流 量 
(m°); 0 为 时 间 步 长 i 内 水 文 站 实际 观测 的 河道 径 
流量 (Gm); QS" 和 0Q 分 别 为 模拟 时 间 段 内 模拟 和 
观测 的 平均 径流 量 (m)。 

2.2 数据 来 源 及 模型 建立 

2.2.1 数据 来 源 与 处 理 驱动 SWAT 模 型 的 数据 有 
DEM 数据 .土壤 数据 .土地 利用 数据 .气象 数据 等 。 
DEM 数据 来 自 地 理 空间 数据 云 平台 ,空间 分 辨 率 为 
90 m(http://www.gscloud.cn/)。 土 壤 数 据 来 自 中 国 科 
学 院 南 京 土壤 所 发 布 的 1:1000000 土壤 数据 库 ,其 
中 土壤 容重 (SOL_BD)、 饱 和 导 水 率 (SOL_K) 和 田间 
持 水 量 (SOL_AWC ) 等 参数 通过 SPAW 模型 计算 得 
到 ,然后 将 数据 重 分 类 成 12 种 土壤 类 型 。 土 地 利用 
数据 来 自 地 理 空间 数据 云 2018 年 30 m 分 辩 率 土地 
利用 数据 (http://www.gscloud.cn/) ,经 过 重 分 类 以 后 
共有 10 种 土地 利用 类 型 。 

水 文 数据 来 源 博 尔 塔 拉 水 文 探 勘 局 ,其 中 包括 
温泉 站 1972 一 2018 年 逐 月 径流 数据 。 气 象 数 据 为 
流域 内 国家 气象 站 点 1968 一 2018 年 的 监测 数据 , 包 
括 降 水 .平均 气温 .最 高 气温 .最 低 气 温 。 用 于 预测 
的 气象 数据 来 自 第 五 次 国际 耦合 模式 比较 计划 
(CMIPS ) 的 预测 数据 ,从 中 选取 8 个 GCM(Global cli- 
mate model) 数 据 ( 表 2) ,包括 中 低 浓 度 排放 情景 
(RCP4.5 ) 和 高 浓度 排放 情景 (RCP8.5 )2 种 情景 模式 
下 2020—2050 年 的 最 高 气温 、 最 低 气 温 、 降 水 数 
据 。 对 这 些 数据 首先 运用 双 线 性 插值 法 将 各 模式 
的 分 辨 率 统一 成 0.5°x0.5° ,然后 采用 RoMBC(Ro- 
bust multivariate bias correction ) 方 法 对 模式 输出 数 
据 进 行 偏差 校正 ,该 方法 能 有 效 降低 数据 的 系统 误 
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表 2 CMIPS 全 球 气候 模式 基本 信息 
Tab. 2 Basic information of CMIPS global climate models 


TRR 二 本 水 平 
编号 。 CMIP5 模 式 。 所 属国 家 WMI phig 
1 ACCESS1-0 澳大利亚 ”CSIRO-BOM 1.9°x1.2° 
2  ACCESSI-3 澳大利亚 ”CSIRO-BOM 1.9°x1.2° 
3 NRM-CM5 法 国 CNRM-CERFACS 1.4°x1.4° 
4 HadGEM2-CC ”英国 MOHC 1.9°x1.2° 
5 HadGEM2-ES ”英国 MOHC 1.9°x1.2° 
6 MIROCS 日 本 MIROC 1.4°x1.4° 
7 MRI-CGCM3 日 本 MRI 1.1°x1.1° 
8 BCC-CSM1-1-M 中 国 NCC 1.1°x1.1° 


Ze) ,最 后 根据 多 模式 等 权重 的 集合 平均 方法 2 输 
出 2 种 情景 模式 下 降水 .气温 数据 ,分 析 2 种 情景 模 
式 下 对 研究 区 径流 变化 的 影响 。 

2.2.2 博 河上 游 SWAT 模 型 构建 本 文 以 0.8 km 作 
为 集 水 国 值 ,将 研究 区 划分 为 19 个 子 流 域 ,1095 个 
水 文 响应 单元 ,在 此 基础 上 输入 历史 气象 数据 构建 
出 SWAT 模 型 。 基 于 水 文 站 实测 月 径流 资料 ,通过 
调整 参数 获得 历史 最 优 模拟 情况 ,运用 模型 评价 指 
标 评 定 模拟 结果 ,最 后 将 CMIP5 多 模式 集合 未 来 情 
景 数据 带 入 率 定好 的 模型 中 ,对 未 来 博 河 上 游 径流 
进行 预测 和 变化 分 析 。 图 2 为 在 构建 SWAT 模 型 模 
拟 方案 时 的 DEM、 子 流域 划分 土地 利用 和 土壤 类 
型 的 空间 分 布 图 。 


(a) DEM 


图 例 


j 
图 例 a Eo 
太 Hk 112 


(c) 土地 利用 类 型 


土地 利用 类 型 人 


故国 耕地 国 国 高 覆盖 度 草地 国 辐 低 覆 盖 度 草地 [水 域 
故国 林地 CO 中 覆盖 度 草地 国 国 永久 性 冰川 雪 地 城镇 用 地 
DANES BAR 
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区 径流 对 气候 变化 的 响应 及 预测 


3 结果 与 分 析 


3.1 模型 参数 化 与 效果 评估 

官方 发 布 的 SWAT 模 型 对 冰川 消融 过 程 模拟 欠 
缺 ,因此 在 模型 构建 时 加 入 冰川 模块 "56。 采 用 swat- 
cup 中 的 SUFI-2 算 法 对 参数 进行 敏感 性 分 析 ,结合 
敏感 性 分 析 结 果 .各 参数 的 物理 意义 以 及 研究 区 特 
殊 的 水 文 过 程 , 采 用 试 错 法 逐步 调整 关键 参数 ,得 
到 最 终 率 定 结果 ( 表 3)。 

以 1968—1971 , 1972—1996 , 1997—2018 年 分 
别 作为 模型 的 预 热 期 . 率 定期 验证 期 ,对 研究 区 月 
径流 进行 模拟 (图 3、 表 4)。 校 准 期 和 验证 期 NSE 分 
别 为 0.81 和 0.83,NSE>0.75, 模 拟 结果 为 优 ;校准 期 
和 验证 期 RSR 分 别 为 0.44 FI 0.41 , RSR<0.5 ,模拟 
结果 为 优 ;校准 期 和 验证 期 PBIAS 分 别 为 0.26% 
和 -6.84% ,IPBIASI<10% ,模拟 结果 为 优 ;校准 期 和 
验证 期 尼 分 别 为 0.82 和 0.86 ,很 接近 1 ,模拟 值 与 实 
测 值 相关 性 很 好 ,模型 整体 评价 为 优 。 以 上 综合 
明 , 加 入 冰川 模块 的 SWAT 模 型 在 博 河 上 游 的 模拟 
结果 较为 可 靠 , 模 型 参数 设置 合理 ,可 用 于 该 流域 
的 径流 过 程 模 拟 和 气候 变化 的 响应 研究 。 
3.2 水 文 气象 因子 变化 分 析 

根据 博 河 上 游 温 泉 站 1968—2018 年 实测 水 文 
资料 及 情景 模式 水 文 资料 ,分 析 其 降水 .气温 变化 


(b) 子 流域 划分 


C | 流域 边界 


(d) 土壤 类 型 


土壤 类 型 ae 
ee ge 石灰 性 黑土 El IT 
ETE meee ee mest 
EE 石灰 性 冲积 士 目 是 永 浆 薄 层 土 国 是 腐殖质 雏形 十 国 松软 潜 育 十 


图 2 数字 高 程 模型 (DEM) 、 子 流域 划分 .土地 利用 类 型 和 土壤 类 型 的 空间 分 布 


Fig.2 Spatial distributions of DEM, subwatershed division, land use types and soil types 
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Tab. 3 Parameter calibration results 
参数 模块 参数 物理 意义 范围 参数 最 优 值 

径流 ALPHA_BF.gw 基 流 a 因子 /d [0, 1] 0.0084 
GW_DELAY.gw 地 下 水 的 时 间 延 迟 /d [0, 500] 120.00 
GWQMIN.gw 发 生 回归 流 所 需 的 浅水 层 的 水 位 阔 值 /am [0, 5000] 100.00 
REVAPMN.gw 渗入 深水 层 所 需 的 含水 层 的 水 位 阔 值 /mm [0, 1000] 1000.00 
CH_K2.rte 主 河道 冲积 物 的 有 效 渗 透 系数 /mm'.h” [0, 500] 30.00 
ESCO.hru LHR AMEE [0, 1] 0.70 

融雪 SFTMP.bsn 降雪 气温 /CC [-5, 5] 1.00 
SMTMP.bsn 融雪 气温 /%C [-5, 5] 0.50 
SMFMX.bsn 6 H 21 日 的 融雪 因子 /mm "qd" [0, 10] 6.50 
SMFMN.bsn 12 H 21 H ÉJRE A/mm Cd" [0, 10] 1.50 
TIMP.bsn KE iat tet I [0.01, 1.00] 1.00 

冰川 Bmelt6 6 H 21 日 的 融 冰 因子 /mm +d [1.4, 16.0] 3.50 
Bmelt12 12 月 21 日 的 融 冰 因子 /mm C!s" [1.4, 16.0] 0.50 
gmlt_tmp 融 冰 温度 阔 值 /Y [-5, 5] 1.44 

3? 一 -实测 径流 一 模拟 径流 

30 上 。 A d 


月 径流 /ms 
a 


1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 2006 2008 2010 2012 2014 2016 2018 
年 份 


图 3 1972 一 2018 年 温泉 站 月 平均 流量 模拟 结果 


Fig.3 Simulation results of monthly average flow of Wenquan station during 1972—2018 


表 4 博 尔 塔 拉 河上 游 流域 模拟 结果 评价 


Tab.4 Evaluation on simulation results of upper watershed of Bortala River 


水 文 站 时 段 NSE RSR PBIAS/% R 
温泉 水 文 站 校准 期 (1972 年 1 月 一 1996 年 12 月 ) 0.81 0.44 0.26 0.82 
验证 期 (1997 年 1 月 一 2018 年 12 月 ) 0.83 0.41 -6.84 0.86 

全 时 段 (1972 年 1 月 一 2018 年 12 月 ) 0.82 0.42 -3.22 0.84 


注 : 尽 为 决定 系数 。 


趋势 (图 4)。2 种 模式 下 平均 气温 分 别 为 5.83 CA 
6.93 C BE BAYH 4.15 % 都 升温 1.50 CHE, Arik 
趋势 显著 ,这 也 与 IPCC 第 一 工作 组 评估 报告 ”预测 
结果 相符 。2 种 气候 模式 下 未 来 降水 量 分 别 为 
328.26 mm 和 302.14 mm, 较 历 史 时 期 年 均 降水 量 
288.63 mm 分 别 增长 了 13.73% (RCP4.5) Fil 4.70% 
(RCP8.5) ,其 增 速 分 别 为 13.21 mm- (10a)! (RCP4.5 ) 


和 4.35 mm: (10a)'(RCP8.5)。RCP4.5 降 水 量 的 增 速 
明显 大 于 RCP8.5 ,其 原因 可 能 是 增 温 导 致 蒸发 量 
加 剧 ,进而 抵消 了 水 循环 加 快 使 降水 增加 的 积极 
影响 ”。 此 外 ,通过 累积 距 平 法 发 现 ,2 种 气候 模式 
的 气温 和 降水 到 2050 年 并 未 出 现 明显 拐点 ,这 表明 
增 温 增 湿 还 将 继续 。 


用 2 种 气候 模式 的 降水 .气温 数据 驱动 率 定 好 
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的 SWAT 模 型 ,分 别 得 到 各 模式 下 的 模拟 径流 。 从 
年 际 变化 趋势 (图 5) 来 看 ,在 气温 和 降水 的 综合 有 影 
响 下 ,RCP4.5 和 RCP8.5 情景 模式 中 ,年 均 径 流量 分 
别 为 4.09x10; m FI 4.36% 10° m’, 较 历史 时 期 年 径流 
Ht 3.16x 10° mw 分别 增加 了 0.93x10 mH 1.2x10° mi， 
增 速 分 别 为 0.31x10* nm.(10aj 和 0.40x10* m?+(10a)' 
冰川 径流 对 总 径流 的 贡献 率 由 历史 时 期 的 27.61% 
提升 到 32.45%(RCP4.5) 和 36.99% (RCP8.5) , BEA 
史 时 期 分 别 增长 了 4.84% 和 9.38% ,冰川 加 速 消融 导 
致 冰川 贡献 率 占 比 提升 是 径流 量 增加 的 主要 原 
。 从 年 内 变化 趋势 (图 6) 来 看 ,RCP4.5 冰川 径流 
主要 是 从 4 月 开始 到 9 月 结束 ,而 RCP8.5 冰川 径流 
主要 是 从 3 月 开始 到 9 月 结束 ,这 也 说 明 升 温 会 导 
致 冰川 消融 时 间 提 前 。2 种 模式 下 模拟 的 月 内 径流 
主要 和 集中 在 7、8 月 ,7 月 达到 峰值 ,峰值 流量 分 别 为 
41.00 m°- s (RCP4.5 ) #1150.65 ms'(RCP8.5), 这 也 
与 冰川 径流 峰值 重合 ,7、8 月 冰川 径流 占 当 月 径流 
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量 65% 以 上 , 占 全 年 冰川 径流 的 60%。 值 得 注意 的 
是 ,冰川 消融 峰值 与 汛期 重合 ,气候 变 暧 导致 冰川 
加 速 消融 会 加 剧 水 系统 脆弱 性 和 河川 径流 的 波动 
性 ,使 汛期 极端 水 文 事件 的 发 生 概率 多 于 历史 时 期 ”。 
3.3 气候 变化 对 径流 的 影响 

将 1972 一 2018 年 实测 数据 和 2020 一 2050 年 2 
种 气候 模式 模拟 数据 的 月 径流 、 月 降水 月 均 温 数 
据 按 照 季节 春 (3 一 5 月 )、 夏 (6 一 8 月 )、 秋 (9 一 11 
A) 冬 (12 月 一 翌年 2 月 ) 划 分 ,分 别 采 用 Pearson 相 
关系 数 进行 双 侧 检验 ( 表 5 、 表 6、 表 7)。 按 照 季节 对 
比分 析 可 以 发 现 , 径 流 在 春季 与 气温 的 相关 性 随 着 
温度 的 上 升 相关 性 也 随 之 上 升 ,这 也 证 实 了 气候 变 
暧 会 导致 冰川 消融 时 间 提 前 ;径流 在 夏季 与 气温 均 
通过 0.01 显著 性 检验 , 且 相 关 性 不 断 提 高 ,这 也 说 
明 夏 季 增 温 会 导致 冰川 消融 加 速 ,使 博 河 上 游 径流 
量 增加 ,这 也 印证 了 迪 丽 努 尔 ' 阿 吉 等 六 得 到 的 结 
论 。 秋 冬 两 季 的 温度 基本 在 0 % 以 下 ,气温 与 径流 


QO O : 
= -e_ 历 史 气 温 从-RCP4.5 气 温 + -RCP8.5 气 温 i) 
ir -o- 历 史 降 水 从 -RCP4.5 降 水 RCP 8. SIE7K 
2 1200 
1968 1973 1978 1983 1988 1993 1998 2003 2008 2013 2018 2023 2028 2033 2038 2043 2048 
年 份 
图 4 温泉 站 水 文 数据 及 情景 模式 年 际 变 化 趋势 
Fig. 4 Inter annual variation trends of hydrological data and scenario models of Wenquan station 
6 = 90 
一 一 实测 径流 一 一 RCP4.5 径 流量 
z 一 -RCP8.5 径 流 一 一 冰川 贡献 率 A n 80 
一 一 RCP4.5 冰 川 贡 献 率 一 -RCP8.5 冰 川 贡献 率 A AA 大 À A 
A br AM & aha ae 170 
Z pe 460 $% 
部 3 E 
2 450 
a = 
迎 2 Jao ® 
oy 
1 i i Á | l pa he +30 
0 20 


1972 1977 1982 1987 1992 1997 2002 2007 2012 2017 2022 2027 2032 2037 2042 2047 
年 份 


图 5 实测 径流 .模拟 年 径流 、 冰 川 页 献 率 年 际 变 化 趋势 


Fig. 3 Inter annual variation trends of measured runoff, simulated annual runoff and glacier contribution rates 
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ERCP4.5 模 拟 径流 = RCPO4.5%k)1]4%3 日 进一步 的 缩减 .; 调节 能 
RCP SHUEY = RCRO4 SIEN 面积 进一步 的 缩减 ,冰川 调节 能 力 下 降 。 
-从 -RCP4.5 冰 川 贡献 率 - * -RCP8.5 冰 川 贡 献 率 
= 90 x, * 
60 a 4 讨论 
507 a 70 Š 
Ta 40b A a 60 %¥ E TE R ah ope 5 a 
$ \ 50 = 博 尔 塔 拉 河上 游 属 于 高 寒山 区 ,径流 主要 以 高 
iw 30 上 会 一 -一 人 al ; = BA a Pelie 
型 ol 山 降 水 和 冰川 融雪 为 主 。 众 多 学 者 在 研究 博 河 径 
w% A ape 196 去 流 时 ,多 是 考虑 降水 和 融雪 对 径流 的 影响 ,而 忽 
视 了 高 寒山 区 冰川 消融 对 径流 的 影响 ,二 现 勇 等 " 


0 
1 2 3 4 5 6 7 8 9 10 11 12 


A} 采用 SWAT+CMADS 模 拟 精 博 河 发 现 温 果 站 的 模拟 
精度 低 于 精 河 站 ,猜测 是 因为 温 果 上 游 高 寒山 区 冰 

图 6 模拟 情景 径流 .情景 冰川 径流 、 冰 用 a AE 
贡献 率 年 内 变化 趋势 川 对 温 果 站 模拟 结果 影响 较 大 。 本 文 将 SWAT 模 型 
Fig.6 Annual change trend of simulated scenario runoff, 和 冰川 模块 耦合 ,该 模块 在 众多 的 以 冰川 径流 占 比 
Scenario glacier runoff and glacier contribution rates 较 高 的 径流 模拟 中 均 有 较 好 的 适用 性 "5223 ,弥补 


了 SWAT 模 型 冰川 消融 过 程 模拟 欠缺 的 问题 ,能 有 
的 相关 性 不 高 ,气温 与 降水 的 相关 性 较 高 ,秋冬 两 。 效 提高 模型 在 高 寒山 区 的 模拟 精度 。 
季 主 要 通过 降水 的 形式 促进 冰川 的 积累 。 但 秋季 气温 变化 是 导致 冰川 消融 速率 改变 的 主要 因 
气温 与 降水 的 相关 性 随 着 气温 的 升 高 ,相关 性 也 随 。” 素 。 通 过 对 未 来 的 气温 分 析 发 现 , 未 来 气温 平均 升 
之 降低 ,这 改变 了 冰川 的 积累 .消融 规律 ,导致 冰川 高 15%。 有 学 者 预 估 了 在 1.5 升温 条 件 下 ,亚洲 


表 5 历史 水 文 数据 突变 对 比分 析 


Tab. 5 Comparison and analysis of abrupt changes in historical hydrological data 


历 由 春季 夏季 
降水 气温 径流 降水 气温 径流 
降水 量 1.000 0.541" 0.013 1.000 -0.143 0.178" 
气温 0.541" 1.000 -0.237” -0.143 1.000 0.503” 
径流 0.013 -0.237 1.000 0.178- 0.503" 1.000 
Ft 秋季 冬季 
降水 气温 径流 降水 气温 径流 
降水 量 1.000 0.360" 0.178" 1.000 0.253” 0.017 
气温 0.360" 1.000 -0.173° 0.253" 1.000 -0.060 
径流 0.178- -0.173 1.000 0.017 -0.060 1.000 


注 :*、** 分 别 表示 通过 0.05 .0.01 显著 性 水 平 检验 。 下 同 。 


表 6 RCP4.S 情 景 模 式 径 流 ` 气 温 、 降 水 相关 系数 


Tab. 6 Correlation coefficients of runoff, temperature and precipitation in RCP4.5 scenario model 


春季 夏季 
RCP4.5 poe = - 一 一 
降水 气温 径流 降水 气温 径流 
降水 1.000 0.704” 0.321" 1.000 -0.004 -0.120 
气温 0.704” 1.000 0.486" -0.004 1.000 0.733" 
径流 0.321” 0.486” 1.000 -0.120 0.733” 1.000 
秋季 冬季 
RCP4.5 - 一 一 - — 
降水 气温 径流 降水 气温 径流 
降水 1.000 0.004 0.103 1.000 0.124 -0.062 
气温 0.004 1.000 0.028 0.124 1.000 0.042 
径流 0.103 0.028 1.000 -0.062 0.042 1.000 
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R7 RCP8.S 情 景 模式 径流 气温、 降水 相关 系数 
Tab. 7 Correlation coefficients of runoff, temperature and precipitation in RCP8.S scenario model 
春季 夏季 
RCP8.5 y Fy >: » 4 y 
降水 气温 径流 降水 气温 径流 
降水 1.000 0.641” 0.236" 1.000 -0.395" -0.155 
气温 0.641” 1.000 0.540” -0.395" 1.000 0.870" 
径流 0.236 0.540" 1.000 -0.155 0.870" 1.000 
RCP8.5 - - 一 一 
降水 气温 径流 降水 气温 径流 
降水 1.000 -0.135 -0.045 1.000 0.261" 0.142 
气温 -0.135 1.000 0.203 0.261° 1.000 0.113 
径流 -0.045 0.203 1.000 0.142 0.113 1.000 


Fey LL RBI 21 PHA AK) R Hae A EO ,这 
意味 着 因 冰 川 消 融 速率 加 快 带 来 的 径流 量 增 加 是 
不 可 持续 的 ,冰川 调节 能 力 的 下 降 会 增加 径流 量变 
RO ,给 水 资源 管理 者 带 来 新 的 挑战 。 

本 研究 仍 存在 不 足 之 处 。 研 究 区 枯水期 相 较 
于 丰 水 期 拟 合 较 差 ,导致 枯水期 模拟 效果 较 差 的 原 
因 是 , 博 河 自 西向 东 流 至 温泉 时 ,经 过 断 陷 盆地 ` 基 
地 隆起 等 地 质 结构 ,导致 地 表 水 地 下 水 频繁 转换 ”9 ， 
频繁 的 地 表 水 和 地 下 水 转换 ,使 得 枯水期 的 径流 量 
呈现 出 无 规则 的 波动 ,进而 增加 模拟 枯水期 径流 量 
时 的 难度 。 今 后 研究 中 考虑 将 SWAT 模 型 与 MOD- 
FLO 允 模型 等 地 下 水 模型 进行 耦合 ,进一步 提高 模 
拟 精 度 。 冰 川 模块 中 的 冰川 消融 因子 ( 冰 辐 射 融 化 
因子 ) 属 于 经 验 参数 ,模拟 径流 时 是 以 多 年 径流 过 
程 线 为 参照 物 进行 调 参 的 ,单纯 的 径流 校正 会 因 降 
水 径流 模拟 的 不 准确 导致 冰川 消融 因子 的 低估 和 
高 佑 ,从 而 影响 径流 组 分 的 分 割 中 。 如 何 降低 因 调 
参 带 来 的 模型 模拟 差异 ,从 而 做 到 精准 化 分 析 径 流 
组 分 是 今后 值得 讨论 的 问题 之 一 。 


5 结论 


本 文 以 博 尔 塔 拉 河 源流 区 为 对 区 ,构建 通信 冰 
川 模 块 的 SWAT 模 型 ,在 CMIP5 不 同情 景 模式 下 , 预 
测 径流 变化 趋势 ,评估 气候 变化 对 径流 的 影响 。 主 
要 结论 如 下 : 

(1) 考虑 冰川 融化 对 径流 的 影响 ,构建 般 入 冰 
川 模块 SWAT 模 型 。 结 果 表 明 , 般 入 冰川 模块 的 
SWAT 模 型 ,全 时 段 的 NSE 为 0.82, RSR 为 0.42， 
RBIAS 为 -3.22% ,RR 为 0.84, 模 型 评价 结果 为 优 ,证 
明 本 次 模型 的 构建 在 该 流域 有 很 好 的 适用 性 。 


(2) 在 CMIP5 的 RCP4.5 和 RCP8.5 情 景 模式 中 ， 
未 来 平均 气温 分 别 为 5.83 % 和 6.93 % ,未 来 降水 量 
分 别 为 328.26 mm 和 302.14 mm, 两 者 较 历史 时 期 均 
表现 出 增加 的 趋势 。 用 2 种 气象 数据 模拟 径流 ,得 
到 的 年 均 径 流量 历史 时 期 年 均 径 流量 分 别 增 加 了 
0.93 10° 呈 和 1.2x10:m。 冰 川 径 流 对 总 径流 的 贡 
献 率 分 别 增长 了 4.84% 和 9.38% ,冰川 加 速 消融 导致 
冰川 径流 占 比 增 大 是 径流 增 大 的 主要 原因 。 
RCP8.5 比 RCP4.5 冰川 消融 提前 了 一 个 月 ,表明 升 
温 会 导致 冰川 消融 时 间 提 前 。 冰 川 消融 峰值 和 汛 
期 重合 ,会 加 剧 了 水 系统 脆弱 性 和 河川 径流 的 波动 
性 ,使 汛期 极端 水 文 事件 的 发 生 概 率 多 于 历史 
时 期 。 

(3) 按照 季节 将 历史 和 2 种 气候 模式 的 径流 、 
降水 .气温 进行 分 类 并 进行 相关 性 分 析 发 现 , 随 着 
气温 的 升 高 ,冰川 消融 时 间 提 前 ,冰川 消融 加 速 , 冰 
川 的 积累 时 间 减 少 ,导致 冰川 面积 进一步 的 缩减 。 
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Response and prediction of runoff to climate change in the 
headwaters of the Bortala River 


TIAN Haowei'’, CHEN Fulong, LONG Aihua™, LIU Jing™, HAI Yang’ 
(1. College of Water and Architectural Engineering, Shihezi University, Shihezi 832000, Xinjiang, China; 2. Department of Water 


Resources, China Institute of Water Resources and Hydropower Research, Beijing 100038, China; 3. State Key Laboratory 


of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China) 


Abstract: Glacial runoff is a major component of runoff in the northwest arid zone of China. Understanding the 
impact of climate change on glacial runoff is crucial, but few studies have been conducted in this field of study in 
the Bortala River Basin, Xinjiang, China. In this paper, we present the glacier module that was added to the 
SWAT model and used to simulate monthly runoff in the headwater area of the upper Bortala River Basin. We suc- 
cessfully simulated monthly runoff at the Wenquan hydrological station during the period 1972—2018. Further, 
we investigated the impact of future climate change scenarios (RCP4.5 and RCP8.5 scenarios for 2020—2050, 
based on CMIP5 climate data) on glacier runoff. The model was able to accurately simulate changes in the source 
area’ s runoff process. The results showed that: For the whole simulation period, the Nash- Sutcliffe efficiency 
was 0.82, the percent bias was —3.22%, the ratio of root mean square error to standard deviation of measured val- 
ue was 0.42, and the coefficient of determination was 0.84, thus allowing the model to be rated as “excellent”. In- 
creasing runoff trends were identified in the simulations of both future climate scenarios, with total runoff increas- 
es of 0.3110* m°» (10a)! and 0.40x10ë m°- (10a)` and increases in the percentage of glacial runoff of 4.84% and 
9.38% , respectively, when compared with the historical period, in which the glacial runoff percentage was 
27.61%. These increases in glacial runoff percentage are the main causes of the increases in runoff volume. Corre- 
lation analysis revealed that as the temperature increases, glacier ablation advances and accelerates, and glacier 
accumulation time decreases, leading to further future shrinking of glacier area. The study provides a basis for 
making changes to historical hydrological information, exploring future evolutionary trends, and mitigating poten- 
tial climate change risks in the region. 


Key words: SWAT model; climate change; CMIP5 climate scenario model; glacial runoff; Bortala River 


